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Abstract 

We explore the abundance of light clusters in asymmetric nuclear matter at subsaturation density. 
With increasing density, binding energies and wave functions are modified due to medium effects. 
The method of Albergo, Costa, Costanzo and Rubbino (ACCR) for determining the temperature 
and free nucleon density of a disassembling hot nuclear source from fragment yields is modified to 
include, in addition to Coulomb effects and flow, also effects of medium modifications of cluster 
properties, which become of importance when the nuclear matter density is above 10~ 3 fm -3 . We 
show how the analysis of cluster yields, to infer temperature and nucleon densities, is modified 
if the shifts in binding energies of in medium clusters are included. While, at low densities, the 
temperature calculated from given yields changes only modestly if medium effects are taken into 
account, larger discrepancies are observed when the nucleon densities are determined from measured 
yields. 

PACS numbers: 25.70.Pq, 21.65.-f, 21.45.-v, 05.70.Ce, 24.10.Pa 
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I. INTRODUCTION 



Understanding of nuclear matter at extreme conditions is one of key issues to clarify 
the problems in core-collapse supernovae as well as neutron stars and heavy ion collisions. 
While heavy ion reactions are often employed to explore the nuclear matter equation of 
state (EoS), careful theoretical work is needed to analyze the experimental signatures and 
to reconstruct the properties of hot and dense matter from the detected abundances and 
energy distributions of ejectiles or from correlation functions between different ejectiles, 
produced in those reactions. 

Among the interesting observables are the yield ratios of different fragments measured in 
such reactions. In many experiments one commonly observes light elements such as neutrons 
(n), protons (p), deuterons (d), tritons (i), 3 He (h), and 4 He (a) (see for example Ref. [l| 
and Refs. therein). Larger clusters, typically with mass numbers 5 < A < 20, are also 
observed and the production process of these fragments must be explained. 

The decay of highly excited nuclear matter produced in heavy ion collisions is a complex 
dynamic process which needs, in principle, a sophisticated treatment. One simple approach 
is the freeze-out concept in which the hot and dense matter in the initial stage is assumed 
to reach thermal equilibrium as long as reaction rates are high. With decreasing density, 
the reaction rates decrease and the equilibration process becomes suppressed. At that time 
the nuclear thermal and chemical equilibrium is frozen out. Often the description of the 
nuclear matter, in particular the distribution of clusters, is calculated within a statistical 
multifragmentation model assuming nuclear statistical equilibrium (NSE) 

as 

. Under the 

simplifying assumption that the final reaction product distribution is identical to the cluster 
distribution at the freeze out point, the thermodynamic parameters such as temperature 
T and particle number densities, n n and n p for neutrons and protons, respectively, can be 
reconstructed from the observed abundances. A simple method for extracting the temper- 
ature of the fragmenting hot system was given by Albergo, Costa, Costanzo and Rubbino 
(ACCR) [4]. The method is based on selecting double isotope (or isotone) ratios, i?2, such 
that the nucleon chemical potentials are eliminated leading to a relation between R 2 , T and 
the binding energies of the isotopes (isotones). This method has been used in the analysis 



of a large number of experiments. See, for example, the early works of Refs. 



labj. 



However, one has to be aware that the dynamic reaction processes do not cease abruptly, 
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so that the concept of a unique freeze out time is only approximate. A single freeze out 
time, independent of the species under consideration and their dynamical state, the flow and 
further parameters describing non-equilibrium effects, may not exist. In addition, secondary 
(post freeze-out) decay will modify the original distribution which also contains excited 
states, and the decay products will be found in the final distribution. We point out that in 
Ref. [81] , the ACCR method was modified to account for the screening due to the Coulomb 
interactions among fragments in the freeze-out volume by using the Wigner-Seitz approxi- 
mation [9]. It was found that the corrections for the temperature are less than 20%, though 
for certain isotone double ratios it can be as large as 50%. In Ref. [10| the ACCR method 
was modified to account for the effect of radial collective flow. It was found that the effect 
on the extracted temperature is relatively small, but the increase in the freeze-out density 
can be significant for large flow energy. It was noted in Ref. |(| that an important improve- 
ment of the simple ACCR method resulted from taking into account post emission decay 
(secondary decay) processes of particles and, in particular, 7 which modify the freeze-out 
yield ratios. Without this correction, different double ratios R2 associated with selected sets 
of fragments (different thermometers) may result in significantly different temperature T. 
We will not discuss these issues further in this work, for a review see for example Ref. 11]. 

If the freeze-out density is not very low, i.e. not at baryon number densities ns < 10~ 4 
fm -3 , the NSE will be modified by medium effects. In this paper, the shift of the binding 
energies of the light elements in hot and dense nuclear matter Q] is considered. The 
changes of the mass fractions of different nuclei due to medium effects complicates the 
determination of the temperature T and baryon number density ns = n n + n p from yields 
of the observed products. In the following we show that simple NSE approaches like the 
ACCR thermometer can be improved if in-medium effects are taken into account. 



II. THE YIELD THERMOMETER 

Within a quantum statistical approach to nuclear matter HQ , using a cluster decom- 
position for the self-energy we obtain expressions for the total proton density 



(T, n P , fi n ) = ^Zn{A,Z) (1) 

A,Z 
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and for the total neutron density 

n n (T,n P ,n n ) = y %2(A-Z)n(A,Z). (2) 

A,Z 



Here, n(A, Z) is the contribution of the A-nucleon cluster to the total nucleon density. Both 
equations <^ and (jSJ) may be considered as a nuclear matter equation of state (EoS) which 
determines the nucleon densities n T as functions of the temperature and the neutron and 
proton chemical potentials, respectively denoted by T, fi n and Additional thermody- 
namic quantities such as free energy and other thermodynamic potentials are obtained by 
integration. 

Starting with the ideal mixture of different species, where the interaction between the 
species is neglected, the number density of a cluster (A, Z) is given by 

n^(A, Z) = g AZ J ^f A> z[E%(p)\, (3) 
with the (Fermi or Bose) distribution function 

Ia,z(E) - ^ [E _ Z/ip _ {A „ z) ^ n] _ TZjya' ^ 

In Eqs. (El), Ef z {j>) = Ef z + H 2 p 2 /(2Am), where, in the non-interactine case considered 
here, E^ z is the ground state binding energy, g^ z is the degeneracy factor of an isolated 
nucleus with mass number A and charge number Z, m is the average nucleon mass and 
H n and fi p are the chemical potentials of neutrons and protons, respectively. Note that in 
general excited states which are characterized by internal quantum numbers may occur, in 
addition to {A, Z}. In that case, a summation over the excited states, including scattering 
states, should be carried out. The number density is assumed to be proportional to the 
cluster yield, observed after freeze-out. In the non-degenerate limit we have the prediction 
of the cluster yields within NSE models. 

Z) oc n^(A, Z) = g AZ ^ e -^- z ^ A ~ z ^' T . (5) 

Inserting n^\A, Z) for the cluster densities n(A, Z), Eqs. (0Q) and ([2]) read 



(T )/V // n ) = ^Zn(°)(A^), 

A,Z 

n\T,fj, p ,n n ) = Y^{A-Z)nW(A,Z), (6) 



A,Z 



for the total proton and neutron densities, respectively. They are approximations to the 
nuclear matter EOS, reflecting NSE. 

Let us now consider the observed cluster yields Y(A, Z) which are proportional to the 
number density fractions n(A, Z)/ub of the cluster {A, Z}. We introduce the (single) ratio 
of the observed cluster yields 

_ Y(A,Z) 

H(AZUA'Z') ~ Y (A',Z'Y ^ 

If we accept the concept of NSE, identifying the observed cluster yields with the predicted 
ones, Eq. (J5j), we can get an estimation for the temperature T and the chemical potentials 
fi n , Hp of nuclear matter produced in heavy ion collisions in a fashion similar to that employ- 



ing the well-known Saha equation in plasma physics 14j . Specifically, since the abundances 
of different bound states are determined by the temperature and the chemical potentials, 
observed yield ratios can be used to determine these parameters. A simple method to derive 
the temperature of the hot system was given by ACCR [4j, assuming NSE (jSJ) and selecting 
double isotope ratios such that the nucleon chemical potentials are eliminated. In particular, 
the H-He thermometer considers the double ratio R^ e of cluster yields Y^°\ 



y(o)(3 H )y(o)(3 He ) 2-2V3-3, 
where the degeneracy and mass factors are explicitly included. Identifying the double ratio 

flHH . = ^™£2L (9) 

-K( 3 He),( 4 He) 

of the observed cluster yields for c?( 2 H), t( 3 H), h ( 3 He) and a ( 4 He) with the prediction 
according to the NSE, -RHHe = -Rhhc' we deduce the ACCR temperature T^ e (= T^ e ) 
corresponding to the observed double ratio -RHHe as 

(a) _ 14.325 MeV 
iHHe -ln[1.591itlHHer [ j 



The constants 14.325 MeV = -[-2.225 - 28.3 + 8.482 + 7.718] MeV and 1.591 = 9/V32 
reflect the ground states binding energies, spins and mass numbers of the ejectiles as given 
in Eq. ©. 

Other combinations of isotopes can be used to construct double ratios of cluster yields 
where, within a simple NSE, the chemical potentials cancel out so that an ACCR temper- 
ature can be derived directly. Thus, thermometers based on the yields of other nuclei such 



as lithium or beryllium isotopes can be introduced. Similar approaches are used in hadron 
production to derive the temperature for the quark-gluon plasma phase transition [l^, or 
in plasma physics [l4, 1(| considering spectral line intensities of different ionization states 
of radiating atoms. 

The advantage of the double ratio is that, within NSE, it does not contain the density, 
because the chemical potentials cancel. Therefore the temperature determination seems to 
be insensitive with respect to the determination of other parameters. These other parame- 
ters, in particular the chemical potentials, are observed if, in addition to the double ratios 
like -Rime, the single ratios R[az),{A'Z') — Y(A, Z)/Y(A', Z') of yields are considered. 

There are some objections to inferring the parameter values of hot dense matter from 
the cluster yields. First, we have to take into account that collisions lead to initially in- 
homogeneous system evolving in time. Even assuming local thermal equilibrium, one has 
to separate the ejectiles arising from different sources. In efforts to do this, the H-He ther- 
mometer has been applied to the double ratio R Vsurt of cluster yields Y(A, Z) for clusters 
with the same surface velocity In that case an additional factor (9/8) arises in the 
temperature equation when the number densities as a function of velocity are employed. 

An important improvement of the simple NSE model was to take into account secondary 
decay processes which modify the freeze-out yield ratios. This has been considered in differ- 
ent papers. This correction is essential to reduce the differences of the ACCR temperatures 
obtained from different thermometers Q]. 

The simple NSE is based on a chemical picture considering a non-interacting, ideal mix- 
ture of different components, which is in chemical equilibrium due to reactive collisions as 
described by the mass action law. Such an approach is valid in the low-density limit, and 
related expressions such as virial expansions can be taken as a benchmark in that limit 



17l | . With increasing density, modifications arise which are based on taking the interactions 



between the different components into account. Thus, as the density increases, corrections 
to the ACCR approach to derive the temperatures of hot and dense matter are expected. In 
earlier work, the effects of the screening of the Coulomb interaction and of the flow on the 



freeze-out density and temperature of disassembling hot nuclei have been considered [lOUllj. 
The effect of screening of the Coulomb interaction becomes of importance for heavy nuclei 
at densities near to the saturation density. 

The main topic we address in this paper is the required modification of the description 
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of the matter as an ideal, noninteracting mixture of different components when densities 
are not low enough to justify this assumption. Despite the fact that the nucleon-nucleon 
interaction is short-ranged, the interaction between the free nucleons as well as nucleons 
bound in clusters is negligible only below about 10 -3 times the nuclear saturation density, 
i.e. at baryonic densities < 10 -4 fm -3 . An important question is the role of medium 
effects due to the nucleon-nucleon interactions. In fact, our work indicates that the concept 
of the simple NSE considering hot and dense nuclear matter as an ideal mixture of different 
clusters is not appropriate to describe disassembling hot matter at densities at and above 
approximately one tenth of saturation density. We address this in the following section. 



III. MEDIUM MODIFICATION OF CLUSTER PROPERTIES 



Recent progress in the description of clusters in low density nuclear matter I2I, Jj], Q, [2^] 
enables us to evaluate the abundance of deuterons, tritons and helium nuclei in a micro- 
scopic approach, taking the influence of the medium into account. Within a quantum statis- 
tical approach to the many-particle system, we determine the single-particle spectral func- 
tion, which allows calculation of the density of the nucleons as a function of T, /i n and [i v . 
The main ingredient is the self-energy E(l,z) which is treated in different approximations. 
The single-particle spectral function contains the single-nucleon quasiparticle contribution, 
E qu (l) = E^ u z (p) or E qn (p), where r denotes isospin (neutron or proton). The quasiparticle 
energy follows from the self-consistent solution of E qu (p) = h 2 p 2 /(2m T ) + ReS[p, E qu (p)]. 

Expressions for the single-nucleon quasiparticle energy E qu (p) can be given by the Skyrme 
mean-field parametrization 2l| or by more sophisticated approaches such as relativistic mean 



field approaches 



22 1 and relativistic Dirac-Brueckner Hartree Fock 



23l | calculations. In the 



effective mass approximation, the single-nucleon quasiparticle dispersion relation reads 

Er(p) = A£f(0) + ^p 2 + 0(p*) , (11) 

where the quasiparticle energies are shifted by AEf E (0), and m* denotes the effective mass 
of neutrons (r = n) or protons (r = p). Both quantities, AE^ E (0) and m*, are functions of 
T, rip and n n characterizing the surrounding matter. Empirical values for the effective mass 
near the saturation density are different from the nucleon mass. In the low-density region 
considered here, the effective mass may be replaced by the free nucleon mass. For calculating 
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the yields, the quasiparticle shift AE^ E (0) can be implemented in a renormalization of the 
corresponding chemical potentials. 

In addition to the £-like quasiparticle contribution, the contribution of the bound and 
scattering states can also be included in the single-nucleon spectral function, by analyzing 
the imaginary part of S(l, z). Within a cluster decomposition, A-nucleon T matrices appear 
in a many-particle approach. These T matrices describe the propagation of the A-nucleon 
cluster in nuclear matter. In this way, bound states contribute to the EoS, n T = n T (T, fi n , /i p ), 



see Refs. |13l . |24|. In the low-density limit, the propagation of the A-nucleon cluster is 
determined by the energy eigenvalues of the corresponding nucleus, and the simple EoS, (CQ) 
and ([2]), results. 

For the nuclei embedded in nuclear matter, an effective wave equation can be derived 
0, Q] • The A-particle wave function and the corresponding eigenvalues follow from solving 
the in-medium Schrodinger equation 

[E*(l) + ■■■ + E*(A) - EZ(pMa, p (1 ...A) 

+ /(*) - fUW(ij, 13) n <W Wi' ■■■A') = o. (12) 

1'.. .A' i<j k^i,j 

This equation contains the effects of the medium in the single-nucleon quasiparticle shifts 
as well as in the Pauli blocking terms. 

The in medium Fermi distribution function /(l) = {exp[E qu (l)/T — fli/T] + 1} _1 con- 
tains the effective chemical potential /2i which is determined by the total proton or neutron 
density, calculated in the quasiparticle approximation, n T = Y^i /(l)^n,r- It describes 
the occupation of the phase space neglecting any correlations in the medium. In the low- 
density and non-degenerate limit (ji T < 0), assuming the effective mass approximation for 
the nucleon quasiparticle dispersion relation, we eliminate jl r using 



exp[ET{p)/T - fi r /T] + 1 ~ 2 \m*T ) 



1 _ nr_ /2?rfi 2 \ 



/r(p) = — Ti^n-TTF :, Wl1 w 7bF e '^- ( 13 ) 



The solution of the in-medium Schrodinger equation ffT2l) can be obtained in the low 
density region by perturbation theory. In particular, the quasiparticle energy of the A- 
nucleon cluster follows as 

EXz(p) = Ef z + ^ + AEf z ( P ) + AE%f( P ) . (14) 



S 



Additional contributions such as the Coulomb shift AE^ ? 1 (p) , which can be evaluated for 
dense matter in the Wigner-Seitz approximation , lll|, |25|, will not be considered here 
since the values of Z are small and the densities are low. The general formalism also allows 
us to describe pairing or quartetting, but this will not be done here. Disregarding the effects 
due to the change of the effective mass, the self-energy contribution to the quasiparticle shift 
is determined by the contribution of the single-nucleon shift 



AEf z (0) — (A — Z)AE™(0) + ZAE™(0) 



SE, 



-.SE, 



(15) 



Inserting the medium-dependent quasiparticle energies in the distribution functions 
/a^IE'a^p)], Eq. (j4]), this contribution to the quasiparticle shift can be included by renor- 
malizing the chemical potentials fi n and fi p . 

The most important effect on the calculation of the yields of light elements comes from 
the Pauli blocking terms in Eq. (1121) in connection with the interaction potential. This 
contribution is restricted only to the bound states so that it may lead to the dissolution of 
the nuclei if the density of nuclear matter increases. The corresponding shift AE A ^ l {p) can 
be evaluated in perturbation theory provided that the interaction potential and the ground 
state wave function are known. After angular averaging, the Pauli blocking shift can be 
approximated as 

\p) « AE^f(0) e~^r . (16) 



AE A ^ Z 



The shift o 
in density 



the binding energy of light clusters at zero total momentum which is of first order 
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191 ] has been calculated recently [12J. Besides neutrons (n) and protons (p), 



light elements deuterons 2 H, {A, Z} = d, tritons 3 H, {A, Z} = t, hellions 3 He, {A, Z} = h, 
and a-particles 4 He, {A, Z} = a have been considered. The interaction potential and the 
nucleonic wave function of the few-nucleon system have been fitted to the binding energies 
and the root mean square (rms) radii of the corresponding nuclei. The following results (in 
MeV, fm) are obtained for the binding energy shifts. 



Pauli 



AE, 



AE Pauli 



A £Pauli 



38384 



0.39402 e 



0.049418 1+- 



3 p 



AE 



Pauli 

a 



(1 + 22.52)1/2 

3389.7 [1 + 0.13347T]~ 3/2 

3901.5 [1 + 0.16455 T]~ 3/2 
4716.0 [1 + 0.09372 T]' 3/2 (n p + n n ) . 



'Erfc 



0.2223 1 + 



22.52 



1/2' 



Tip H~ Th n 

T s / 2 



n, 



4 2 
-n„ -I — n T 
3 p 3 



(17) 



These results describe only the linear shifts as functions of the nucleon densities. The 
differences between the values for AEf auh and AE^ anh are mainly caused by different values 
of the rms radii for these two nuclei. With increasing density, higher orders terms with 
respect to the densities also become relevant. 

It can be shown [25] that the EoS can be evaluated as in the non-interacting case given 
above, except that the number densities of clusters must be calculated with the quasiparticle 
energies, 

n^(A, Z) = g AZ J ^- 3 fA,z[EXziP)] • (18) 
In the cluster-quasiparticle approximation, the EoS, (0Q) and (j2J), reads 



A,Z 

< U (T, fi p , //„) = J2( A - Z)n qn (A, Z) , (19) 

A,Z 

for the total proton and neutron density, respectively. 

This result is an improvement of the NSE and allows for the smooth transition from the 
low-density limit up to the region of saturation density. The bound state contributions to 
the EoS fade with increasing density because they merge with the continuum of scattering 
states. This improved NSE, however, does not contain the contribution of scattering states, 
in particular resonances appearing in the continuum of scattering states when bound states 
merge with the continuum. For the treatment of scattering states in the two-nucleon case, 



as well as the evaluation of the second virial coefficient, see Refs. |17l . I24j . We will also not 
consider the formation of heavy elements here. This limits the present results to the range 
of parameters T, n n and n p , where the EoS is determined only by the light elements. For a 
more general approach to the EOS which takes also the contribution of heavier clusters into 



account see Ref. [25]. 



IV. IMPROVED THERMOMETER AND DENSITY DETERMINATION IN- 
CLUDING MEDIUM EFFECTS 

To show the effect of in-medium corrections, we start with a temperature T and densities 
n n and n p and calculate the corresponding yields Y qu {A^ Z), taking the in-medium shifts into 
account. Then we use these yields to infer the parameter values T^ a \= T^), rvh (— n n ) and 
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n p a \= n^), using the ACCR relations which were derived neglecting in-medium corrections. 
In this way we obtain for given ratios of cluster yields R(a,z),{A',z'), Eq. ([7]), {T,n n ,n p } 
which we identify as the values which would be derived from experiments if the medium 
effects are considered and those {T^, ni?\ n p a ^}, derived without taking the medium effects 
into account. Comparing these sets of the parameters we demonstrate how the medium 
modification of the binding energy of light nuclei, Eq. (1181) can modify the results determined 
from the experimental yields of light clusters. The three ratios necessary to determine 
three thermodynamic parameters are derived here from the four yields of the light clusters 
F( 2 H), F( 3 H), F( 3 He) and F( 4 He). 

We first compare results of the determination of the temperature T from cluster yields, if 
the in-medium quasiparticle shifts are taken into account, with the temperature deter- 
mined from the same yields if medium effects are neglected. In particular, the temperature 
T H He is not related in a simple way to the double ratio -RHHe, Eq- ©• Considering the yields 
in quasiparticle approximation, we have 

y q u (2H)yq u ( 4 He) 

yqu(3 H )y q u(3 He )- ^ v ) 

If we take the yields Y qu (A, Z) oc n qu (A, Z), Eq. (|T8l) . in the non-degenerate case, we obtain 
the relation 

T hhc = : M * „ 7 (14.325 MeV + AE^ uli + AE^ uli - AEf auli - AE^ uli ) , (21) 

ln[1.591i?HHeJ 

The energy shifts are functions of temperature and densities so that this relation has to be 
solved self-consistently. 

On the other hand, neglecting in-medium corrections, we find from the same double ratio 
the apparent ACCR temperature T^ e according to Eq. ffTUj) . Using, in the non-degenerate 
case, Eq. ([5]), the relation between both quantities is given by 

-i 



AEf auli + AE^ uh - AE% auli - AE^ uh 

1 H P (0) , p (0) ^(0) ^(0) 

E t + E h - E d - E a 



TuHe ■ (22) 



In the approximations considered here The self-energy contributions to the shifts disappear, 
in a fashion similar to the chemical potentials. In Figure 1 we show the ratio between the 
ACCR temperature and T as a function of the baryon density n# = n p + n n for various 
values of T. 
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Similarly, the densities can be estimated by considering single ratios R(a,z),(A',z'), Eq. 
(J7|). If the shifts of the binding energies due to medium effects are neglected, we have 

9A,zA V2 



R 



(a) 

{A,Z),{A',Z>) 



7 (0) _ 



E^> ,-(Z-Z')^ a) -(A-Z-A'+Z')^ n a) ]/T^ 



(23) 



9A,z'A'^ 

Assuming NSE and considering special combinations, we can obtain the chemical potentials 
of protons (fj, p ) and neutrons (/i n ) from the triton to 4 He ratio or from the 3 He to 4 He ratio, 
respectively, as 



(a) 



-19.8MeV + T (a) In 



-20.6MeV + T (a) In 




(24) 



(25) 



_ 2 2 Y h _ 

This allows us to calculate the chemical potentials separately. Then, considering chemical 
equilibrium between the different clusters, the total proton and neutron densities are given by 
the mass action law, cf. Eqs. (TTJ) and ([2]). Assuming NSE where in-medium corrections are 
neglected, we find from Eq. ([6]) the total densities n p a \T^ a \ /4 , fj^/) and nl"' '(T^, /4> , /i„ ) 
of protons and neutrons, respectively. 

Taking into account the in-medium quasiparticle energy shifts of the nuclei, the relations 
are changed so that 

'3 3 / 2 Y a - 



a = -19.818MeV + AE* U - AE? U + Tin 



H n = -20.582MeV + AE% U - AE% U + Tin 



2 2 Y t 

3 3 / 2 Y n 



(26) 



(27) 



2 2 Y h _ 

where the temperature T is obtained taking the medium modifications of the energies of 
nuclei into account. Now, the total proton and neutron densities are calculated from the 
EoS, Eqs. ( Tl9|) . which contain medium-dependent quasicluster energy shifts. 

To show the effect of these medium modifications, we start with given values for T, n p and 
n n and calculate the cluster abundances solving Eqs. (fT9l) . taking the shifts into account and 
restricting our consideration to A < 4. This gives us certain values for the chemical potentials 
Hp and fi n . Obviously, within a self-consistent calculation we can not only reproduce the 
input quantities n p and n n from these values of \x v and //„, but also the single ratios for 
different yields, in particular i?4 He 3 H , -R4 He 3 H e, and the double ratio i?HHe- Now, we consider 
this as input and determine within the simple NSE the ACCR values T^ a \ and (jffl. 
In NSE, where medium shifts are neglected, we calculate the number densities (A, Z) 
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of the nuclei, using the EoS, Eq. (Q. Obviously, the single and double ratios given above 
are reproduced. However, not only the temperature will differ from the input value T, 
but also the total proton density n p a ^ and the total neutron density nffl will deviate from 
the input values n p and n n , respectively. In Figure 2 we show the ratio between the ACCR 
baryon density rig = n p + Tin and % = n p + n n as a function of % for various values of 
T. 

V. DISCUSSION AND CONCLUSIONS 

The assumption of NSE provides a simple means to estimate the thermodynamic param- 
eters of nuclear matter at freeze-out from the observed yields of nuclei. This approach is 
applicable as long as the interaction between the clusters can be neglected. However, the 
thermometers and chemical potentials are no longer correctly scaled when the shifts of the 
binding energies due to the interaction with the surrounding matter become of relevance. 
The derivation of the thermodynamic parameters from the measured yields has to be carried 
out in a self-consistent manner since the binding energies, which determine the yields, are 
themselves dependent on the temperatures and densities. 

Analyzing empirical data, the use of the ACCR method can only give a first approxi- 
mation to the temperature and the density. Taking these first estimations, the shift of the 
binding energies of the clusters can be estimated. With these modified energies, the next 
iteration deriving the values of the parameters from the measured yields can be made, and 
a self-consistent solution is expected after a sufficient number of iterations. Alternatively 
one can also produce tables for yields taking the medium shifts into account, so that the 
optimal values of the parameters are obtained by interpolating within the table to identify 
the values of the parameters which best corresponds to the measured yields. 

Comparing the values of the parameters obtained in the full calculation, with inclusion of 
medium effects on the yields with those deduced in the ACCR approach, we find that moder- 
ate deviations in the temperature arise for densities larger than 0.0001 fm -3 . Determination 
of the densities is more sensitive to the medium effects. 

The shift of the binding energies has been given in first order of the density, and higher 
orders terms in the density are expected to contribute if the density increases. Starting with 
baryonic densities near 10~ 2 fm -3 , the composition has to be calculated with momentum 
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dependent shifts instead of the rigid shifts considered here, and then the temperature is 
found from the ratio of the mass fractions after the composition is calculated in a self- 
consistent way. One has to perform the full momentum integration instead of considering a 
rigid shift as given at P = 0, when the shifts depend on the center-of-mass momentum of 
the cluster. The results given here are applicable at densities which are not too high, i.e. 
up to 0.01 fm~ 3 . 

In conclusion we point out that the fragment yields from hot and dense nuclear matter 
produced in heavy ion collisions can be used to infer temperatures and proton/neutron 
densities of the early stages of the expanding hot matter. The assumption of thermal 
equilibrium can be only a first approach to this non-equilibrium process. To determine the 
yield of the different clusters, a simple statistical model neglecting all medium effects, i.e., 
treating it as an ideal mixture of non-interacting nuclei, is not applicable when the density 
is larger than 0.0001 fm~ 3 . Self-energy and Pauli blocking will lead to energy shifts, which 
have to be taken into account to reconstruct the thermodynamic parameters from measured 
yields. The success of the simple ACCR method to derive the values for the temperature 
can be understood from a partial compensation of the effect of the energy shifts so that 
reasonable values for the temperature are obtained also at relatively high densities. More 
care must be taken in inferring densities from the data. It should be mentioned that similar 
questions have to be considered when hadron production is investigated at the quark-gluon 
phase transition. 

Cross-checks can be performed to see to what extent the approach given here is consis- 
tent. Hitherto we considered only the yields of d, t, h and a, and the corresponding ratios 
are reflected by the temperature and the chemical potentials of the neutrons and protons. 
The determination of the yields of additional clusters will allow for a comparison between 
predictions and experimental data. 

Acknowledgments 

S. Shlomo would like to thank the Institut fur Physik of the University of Rostock, 
Rostock, Germany, for the kind hospitality. This work was supported by the US Department 
of Energy under grant No. DE-FG03-93ER40773 and the Robert A. Welch Foundation under 



14 



grant No. A0330. 



[1] S. Kowalski et al, Phys. Rev. C 75, 014601 (2007). 

[2] J. P. Bondorf, A. S. Botvina, A. S. Iljinov, I. N. Mishustin and K. Sneppen, Phys. Rep. 257, 
133 (1995). 

[3] D. H. Gross, Rep. Prog. Phys. 53, 167 (1990). 

[4] S. Albergo, S. Costa, E. Costanzo and A. Rubbino, Nuovo Cimento A 89, 1 (1985). 

[5] J. Pochodzalla et al, Phys. Rev. Lett. 75, 1040 (1995). 

[6] A. Kolomiets et al, Phys. Rev. C 54, R472 (1996). 

[7] J. A. Hauger et al, Phys. Rev. Lett. 77, 235 (1996). 

[8] A. Kolomiets, V. M. Kolomietz and S. Shlomo, Phys. Rev. C 55, 1376 (1997). 

[9] E. Wigner and F. Seitz, Phys. Rev. 46, 509 (1934). 

[10] S. Shlomo, J. N. De and A. Kolomiets, Phys. Rev. C 55, R2155 (1997). 

[11] S. Shlomo and V. M. Kolomietz, Rep. Prog. Phys. 68, 1 (2005). 

[12] G. Ropke, Phys. Rev. C, in production; arXive nucl-th/0810.4645. 

[13] G. R6pke, M. Schmidt, L. Miinchow and H. Schulz, Nucl. Phys. A399, 587 (1983). 

[14] M. N. Saha, Phil. Mag. 40, 472 (1920); Zeit. fur Physik 6, 40 (1921). 

[15] P. Braun-Munzinger, K. Redlich and J. Stachel, in Quark Gluon Plasma 3, Eds. R. Hwa, 

X.-N. Wang, (World Scientific, Singapore, 2004) P. 491 [arXive [nucl-th/0304013"] . 

[16] X. Chen and P. Han, J. Phys. D: Appl. Phys. 32, 1771 (1999). 

[17] C. J. Horowitz and A. Schwenk, Nucl. Phys. A 776,55 (2006). 

[18] G. Ropke, A. Grigo, K. Sumiyoshi and Hong Shen, Part, and Nucl. Lett. 2, 275 (2005). 

[19] G. Ropke, A. Grigo, K. Sumiyoshi and Hong Shen, in Superdense QCD Matter and Compact 

Stars, Ed.: D. Blaschke and A. Sedrakian (Springer, Dordrecht 2006), p. 75. 

[20] K. Sumiyoshi and G. Ropke, Phys. Rev. C 77, 055804 (2008). 

[21] D. Vautherin and D. M. Brink, Phys. Lett. B 32, 149 (1970); Phys. Rev. C5, 626 (1972). 

[22] S. Typel, Phys. Rev. C71, 064301 (2005). 

[23] J. Margueron, E. van Dalen and C. Fuchs, Phys. Rev. C 76, 034309 (2007). 

[24] M. Schmidt, G. Ropke and H. Schulz, Ann. Phys. (N.Y.) 202, 57 (1990). 

[25] G. Ropke, M. Schmidt and H. Schulz, Nucl. Phys. A424, 594 (1984). 

15 



1.2 



1.15 



^ 1.1 h 



1.05 



1 



~i 1 1 1 — i — i — r 



~i 1 1 1 — i — i — r 



0.0001 





.*/ 


■ ■ ■ T = 10 MeV 


T = 8 MeV 




- ■ T = 6 MeV 


* • — 


T = 4MeV 






.*/ /_ 




7 1 




:/ / - 



:/ 

. . 4 

.V / 




J I I I I I L 



0.001 

n fi [fm ] 



0.01 



FIG. 1: (color online) The ratio between the ACCR temperature (no medium effects) and T 
(including medium effects) as a function of the baryon density ub = n p + n n for various values of 
T. 
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